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Parametric solutions of Pell equations 


Leonardo Zapponi 


Abstract. — This short paper is concerned with polynomial Pell equations 

P 2 - DQ 2 = 1, 

with P,Q,D£ C[X] and deg(D) = 2. The main result shows that the polynomials P and Q are 
closely related to Chebyshev polynomials. We then investigate the existence of such polynomials 
in Z[X] specializing to fixed solutions of ordinary Pell equations over the integers. 


Introduction 

This paper is a digestet, condensed and updated version of a discussion recently held in 
the mathematical forum Mathoverflow [1], The original and motivating question, asked by 
Stefan Kohl, can be summarized as follows: 

Let n be a positive integer which is not a perfect square and consider a fundamental 
solution ( a , b ) of the Pell equation 

x — ny =1. 

Is it always possible to construct polynomials P,Q,D £ Z[X], with deg(D) = 2, 
and an integer k, such that 

P 2 - DQ 2 = 1, 

with P(k ) = a , Q(k ) = b and D(k ) = n? In case of positive answer, is it possible to 
bound the degree of P? 

Numerical evidences suggested that given a fundamental solution of the Pell equation, 
the polynomials P , Q and D can always be constructed and that the degree of P is at most 
6. David Speyer posted an explicit family with P of degree 1 and the author did the same 
for P of degree 2, leading to a positive answer to the first question. In a second post, the 
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author was able to prove that, if the integer n is square-free and not congruent to 1 modulo 4 
then the degree of P is at most 2. The general case which was proved subsequently and just 
announced on the forum confirms the numerical evidences as soon as we suppose that n is 
square-free. 

The paper is organized as follows: in the first section, we study the general polynomial 
Pell equation 

P 2 - DQ 2 = 1, 

with D of degree 2. The key result of the paper is theorem 1, which states that P, Q and D can 
be explicitly computed and that they are closely related to Chebyshev polynomials. One of 
the main ingredients of the proof comes from Grothendieck’s theory of dessins d’enfants. The 
first two paragraphs of the second section describe all the families of polynomials P, Q and 
D satisfying the conditions in Stefan Kohl’s question, provided that the degree of P is less 
than or equal to 2. More precisely, for P of degree 1, proposition 2 provides essentially the 
same family obtained by David Speyer, while for P of degree 2, proposition 3 gives a refined 
version of the family posted by the author. The last paragraph is concerned with the possible 
degrees of P. The main result is theorem 5 which asserts that for a fixed solution (a, b), 
the degree is bounded by an explicit constant depending on (a, b ) and corollary 6 establishes 
that without any restriction on the integer n, there is no uniform bound. Nevertheless, if n is 
square-free, corollary 8 claims that the degree is bounded by 6. 
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1. Chebyshev polynomials 

Let K = C(X) denote the field of rational functions over the complex numbers. The set 
of sequences ( u n ) of elements of K satisfying the relation 

tln +1 = 2Xu n U n — \ 

is a A'-vector space of dimension 2, spanned by the sequences (T n ) and (Un ) defined by 

{ To = 1 ’ and 

\Ti=X, \U 1 =2X. 

For any non-negative integer n, the polynomials T n ,U n £ Z[A'] are respectively called 
Chebyshev polynomials of the first and second kind of degree n. It is easily checked that in 
the field L = AT(V X 2 — 1), they satisfy the identity 

(X + VX 2 -l) n = T n + Un-WX 2 - 1, 
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which can be considered as an alternative definition. In particular, the Pell relation 

(1) T 2 - (X 2 - 1 )C/ 2 _ X = 1 

holds. This last identity will be the main ingredient in the following sections. 

From a geometric point of view, the polynomial T n induces a cover of the projective line 
unramified outside oo, 1 and — 1. It follows that the morphism f n : P 1 —> P 1 associated to 
the polynomial 

f=\(T n + 1) 

is a Belyi map. Its ramification data of depends on the parity od n. In any case, it is totally 
ramified above oo. If n is even then all the points above 1 have ramification index 2 and all 
the points above 1 have ramification index 2 expcepted two of them, which are unramified. 
If n is odd then the ramification behaviour above 0 and 1 is the same: all the points have 
ramification index 2 excepted one of them, which is unramified. Grothendieck’s theory of 
dessins d’enfants provides an elegant combinatorial description of the isomorphism classes 
of Belyi maps. We will not enter into the details, but it follows that given a fixed positive 
integer n, there is a unique dessin d’enfant corresponding to the above ramification data (see 
the picture below). We refer to the book [2] for a complete introduction to the subject. 


» 


(n)-•—o 


FIGURE 1. The dessin d’enfant associated to T n . 


2. Polynomial Pell equations 

We now describe the solutions of Pell equation 

P 2 - DQ 2 = 1, 

with P,Q,D £ C[X], in the special case where D has degree 2. The following result shows 
that the polynomials P and Q are closely related to Chebyshev polynomials. 

Theorem 1. — Let P,Q,D £ C[A'] with deg (D) = 2 and deg (P) = d. The following 
conditions are equivalent: 

1. We have the identity 

P 2 - DQ 2 = 1. 

2. There exist constants A, p £ C x and is £ C such that 

(P = ±T d (XX + is), 

< Q = pU d -\(\X + is), 

[D = p~ 2 ((XX + is) 2 - 1) . 
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Proof. — First of all, the relation 1 shows that the second condition implies the first. Sup¬ 
pose now that P,Q,D £ C[X] satisfy condition 1, with deg(D) = 2. Assuming P 0, it 
follows that deg(P) = d > 1 and thus deg(Q) = d — 1. Consider the polynomial / = P 2 , 
so that f has degree 2d — 1. Setting 

r s 

P = - x f) ei and Q = v]J(X - y t ) fi , 

i= 1 2=1 

with u,i)GC,r<d and s < d — 1, we obtain the factorization 

S' = IK* - ^) 2ei_1 t[(X - yi) 2 f '-'R. 

i= 1 i=l 

with R G C[X]. Since d = Y^i =i = 1 + 2i=i /»» we find the identity 

r s 

2d - 1 = ^(2e* - 1) 4- ^(2/i - 1) + deg(i?) = 4d-2-r-s + deg(il), 

2=1 2=1 

which leads to 

r + s = 2d—1 + deg(P). 

It then follows that r = d, s = d—1 and deg(P) = 0, i.e. P and Q are separable. Remark 
that the polynomial D is then itself separable. In this case, the cover P 1 —Y P 1 induced by / 
is only ramified above oo, 0 and 1, i.e. it is a Belyi map. Moreover, it is totally ramified above 
oo, all the points above 0 have ramification index 2 and the points above 1 have ramification 2 
excepted two of them, which are unramified, corresponding to the roots of D. It then follows 
from the discussion in the first section that there exist constants A £ C x and u £ C, such that 

S = \ (r 2d ( AX + V) +1). 

Now, the general identity 

T 2d = 2T 2 - 1 

leads to the relation 

p 2 = t |( xx+ u), 

from which we deduce the expression of P in the theorem. Finally, the polynomials Q 
and D are obtained from relation 1 combined with the fact that the polynomial U d -\ is 
separable. □ 

Remark. — If d is odd, we can remove the ± sign in the expression of P since in this case 
the Chebyshev polynomial T d defines an odd function, so that T d (—XX—i/) = —T d (\X+v). 

3. Parametric solutions 

We are now ready to study Stefan Kohl’s problem stated in the introduction, from which 
we keep the notation. We will refer to the polynomials P,Q,R £ C[X] as a paramet¬ 
ric solution associated to ( a,b ), its degree being the integer deg(P). Remark that we can 
always reduce to the case k = 0. From now on, we restrict to this situation. With this as¬ 
sumption, if P,Q,D £ Z[X\ is a parametric solution, the same holds for the polynomials 
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P(mX), Q(mX) and D{mX) for any integer m ^ 0. In the following paragraphs, we ex¬ 
plicitly answer the first question by performing all the parametric solutions of degree 1 or 
2 having inregral coefficients. We then bound the degree of general parametric solutions, 
showing, among the others, that if the integer n is square-free then the degree is at most 6. 

3.1. Degree 1. — We start by giving an explicit description of all parametric solutions with 
P of degree 1. 

Proposition 2. — Let n be a positive integer which is not a perfect square and consider a 
solution (a, b) of the Pell equation 

2 2 i 

x — ny =1. 

Given three polynomials P,Q,D £ Z[X] with deg(P) = 1 and deg (D) = 2, the following 
conditions are equivalent: 

1. We have the identity 

P 2 - DQ 2 = 1, 

with -P(O) = a, (5(0) = b and D( 0) = n. 

2. Set c = 1 if b is odd and c = 2 if b is even. There exists an integer m 0 such that 
P = Po(c~ 1 mX),Q = band D = Do(c~~ 1 mX), with 

f P 0 = b 2 X + a, 

| P 0 = b 2 X 2 + 2aX + n. 

Proof. — First of all, if the polynomials P, Q and D are defined as in condition 2, it is easily 
checked that they satisfy condition 1. Suppose now that P, Q, D £ Z[A'] satisfy condition 1. 
The polynomial Q is constant, so that Q = b. Following the notation in theorem 1, we obtain 
the identity p = b and, taking account of the remark after theorem 1, we have the expression 

P = XX + v, 

which implies that A is an integer and that v = a. Since b~ 2 {a 2 — 1) = n, we have the 
identities 

D = b~ 2 ((AX + v) 2 - 1) = X 2 b~ 2 X 2 + 2A ab~ 2 X + n, 

and thus b divides A and b 2 divides 2Aa. The integers a and b being coprime, we deduce that 
b 2 divides 2A. Setting A = bu , it follows that b divides cm, say u = c~ 1 bmn, which leads to 
the desired expressions. □ 

3.2. Degree 2. — We now give a complete description of parametric solutions when the 
polynomial P has degree 2. 

Proposition 3. — Let n be a positive integer which is not a perfect square and consider a 
solution (a, b) of the Pell equation 

x — ny =1. 

Given three polynomials P,Q,D £ Z[X] with deg(P) = deg(P) = 2, the following condi¬ 
tions are equivalent: 
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1. We have the identity 

P 2 - DQ 2 = 1, 

with P( 0) = a, Q(0) = b and D(0) = n. 

2. There exist two integers m 0 and e £ {±1} such that, setting 

c = gcd (b 3 } (a + s)b , 2(a + e) 2 ) , 

we have P = Po(c _1 mX), Q = Qo(c~ 1 mX) and D = Do(c~ 1 m.X), with 

(Po = 6 4 (a + e)X 2 + 2 6 2 (a + e)X + a, 

\ Qo = b 3 X + b, 

[l?o = (« + e) 2 b 2 X 2 + 2 (a + e) 2 X + n. 


Proof. — First of all, a direct computation shows that if P, Q and I) are defined as in con¬ 
dition 2, they satisfy condition 1. Suppose now that P, Q and D satisfy the first condition. 
Following theorem 1, we have the expressions 


( eP = 2A 2 X 2 + WvX + 2v 2 - 1, 

\ Q = 2A f-tX + 

[p = A 2 fi~ 2 X 2 + 2A vfjT 2 X + (y 2 — \)y~ 2 . 


with e = ±1. The relations P(0) = a and Q( 0) = b lead to the identities 

ea +1 . b 


v 2 = 


In particular, we obtain the relation 


- and u = —. 

2 2v 


( v 2 -l)y ~ 2 = 


The leading coefficient of Q being an integer, there exists u £ Z such that 



We then deduce that D belongs to Z[X] if and only if b 2 divides m(ea + 1) and b 3 di¬ 
vides 2 u(ea + l) 2 . These last two conditions can be summarized as the divisibility of 
u(ea + 1) gcd(6,2 (ea + 1)) by b 3 , which is itself equivalent to the divisibility of u by c -1 6 3 . 
Setting u = c~ 1 b 3 m, we obtain the expressions in condition 2. □ 


Example. — Setting m = ±1 in the above proposition leads to parametric solutions whose 
coefficients have the smallest absolute value. For example, for n = 31, a = 1520 and 
b = 273, setting m = e = 1, we find c = 59319 and obtain the polynomials 

! P = 2401X 2 + 3822X + 1520, 

Q = 343X + 273, 

D = 49X 2 + 78X + 31, 

whose coefficients are much smaller than those of Po, Qo an d Po, since we have the expres¬ 
sions 

( P 0 = 8448503770161X 2 + 226717218A' + 1520, 

< Qo = 20346417X + 273, 

[p> 0 = 172418444289A 2 + 4626882X + 31. 
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3.3. Bounding the degree. — We finally study the possible degrees of a parametric solu¬ 
tion. We start with an easy but usefull lemma. 

Lemma 4. — Let P,Q, D £ Q[X] be a paramteric solution of degree d associated to (a, b). 
Then there exist Pi, Q\, D\ £ Z[X\ defining a parametric solution of degree d associated to 

(a, b). 

Proof — By assumption, the constant terms of P, Q and D are integers. If m denotes the 
lowest common multiple of the denominators of the coefficients of P, Q and D then the 
polynomials P\ = P(mX), Q\ = Q[mX) and D\ = D{mX) belong to Z[X\ and define a 
parametric solution of degree d associated to (a, b). □ 

Let now n be a positive integer which is not a square and consider the real quadratic 
field K = Q(y/n). Denote by Ok the ring of integers of K. The subgroup U of O^ 
consisting of elements having norm 1 is isomorphic to Z x Z/2Z. The elements of the group 
V = U D Z[y/nJ bijectively correspond to the solutions of the Pell equation 

x — ny =1. 

Moreover, U and V are isomorphic to Z x Z/2Z, the quotient U/V is cyclic and V is gen¬ 
erated by —1 and a + by/n, where (a, b) is a fundamental solution. More generally, if (a, b) 
is any solution of the above Pell equation, denote by V(a, b) the subgroup of V generated by 
a + by/n and —1. 

Theorem 5. — Let n be a positive integer which is not a square and consider a solution 
(a, b) of the Pell equation 

2 2 -i 

x — ny = 1. 

The following conditions are equivalent: 

1. There exists a parametric solution P,Q, D £ Z[X] of degree d associated to ( a, b). 

2. The integer d divides 2 (U : V (a, b)). 

Proof — We start by supposing that the first condition is fulfilled. Let P,Q,D £ Z[X] 
be a parametric solution of degree d. We know from theorem 1 that there exist constants 
A, p £ C x and v £ C such that 

(P = ±T d (\X + v), 

< Q = pU d -i{\X + v), 

[D = p~ 2 ((XX + v) 2 - 1) ■ 

We can restrict to the case where P = Td{XX + v). Indeed, for P = —Td(XX + v), if 
we replace P with — P, we obtain P(0) = —a and Q(0) = 6, which still defines a solution 
of the same Pell equation and the groups V(a, b) and V(—a, b), which are isomorphic, have 
the same index in U. Considering the quotient of the coefficient of X in D by its leading 
coefficient, it follows that A^ _1 is rational. Similarly, comparing the leading coefficients of 
P and Q, we deduce that A/i -1 is also rational. From the expression of the coefficient of X 
in D, the same holds for i '/i -1 . Finally, taking the quotient of the constant term of D by its 
leading coefficient implies that A 2 is rational. If d is odd then the expression of the leading 
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coefficient of P leads to the relation X d £ Q and thus A, v and // are rational. We then have 
the identity 

p + Q\/d = (ax + u + . 

Evaluating at 0, we obtain the relation 

a + b\/n = (y + fi\/n) d , 

In particular, the element v + fiyfn £ K is an algebraic integer; it therefore belongs to U. 
The group U/V(a,b) being cyclic, we deduce that d divides ( U : V{a,b)). Suppose now 
that d = 2 m is even. Setting 


we have the identity 


P 


\P 1 = 4(AX + v) 2 — 1, 
= 2ApA -f- 2r , ^x, 

+ QVD = ^ / ? i + Qi\/1) 


Remark that the above discussion implies that Pi and Q\ belong to Q 1 [A - ]. Proceeding as 
before, we deduce that m divides (U : V(a, b)). 

Let now d be an integer dividing 2 (U : V (a, b)). As before, we start with the case where d 
is odd, so that it divides (U : V(a, b)). From the definition of V (a, 6), there exists an element 
u + v^fn of U such that 

a + b\/n = (u + vy/n) d . 

Remark that 2 u and 2m; are integers (actually, setting n = r 2 s, with s square-free, it follows 
that 2 rv is an integer). In particular, the elements u and v are rational. In this case, the 
polynomials 

'P = T d (X + u), 

Q = vU d _i(X + u), 

D = v~ 2 X 2 + 2 v~ 2 uX + n, 
belong to QfA - ] and define a parametric solution of degree d associated to (a, b). We then 
apply lemma 4. Suppose now that d = 2m is even, so that m divides (U : V(a,b)). As 
before, there exists u + Vyfn £ U such that 

a + b\fn = (tt + v\fn) m . 

Remark that the polynomials P 0 , Q 0 and D 0 defined in roposition 3 also work for rational 
solutions of the Pell equation 

2 2 i 

x — ny =1. 

In the present case, the polynomials 

' Pi = v A (u + 1)X 2 + 2v 2 {u + 1)X + u, 

Qi = v 3 X + v, 

D = (u+ l) 2 v 2 X 2 + 2(u + 1) 2 X + n , 
belong to Q[A] and define a parametric solution of degree 2 associated to (u. v). Setting 

P + Q^D= (Pi+QiVDy , 
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the polynomials P,Q,D £ Q[X] are a parametric solution associated to (a, b) and we apply 
once again lemma 4. □ 


This last result implies that the degree of a parametric solution associated to a fixed solu¬ 
tion (a, b) is bounded. It is then natural investigate the existence of an uniform bound. From 
the above proof, it easily follows that the answer is negative if we don’t restrict to fundamen¬ 
tal solutions. The following result shows that without any assumtion on the integer n, such 
an uniform bound does not exist even when restricting to fundamental solutions. 

Corollary 6. — For any positive integer d, there exist an integer n which is not a square, a 
fundamental solution of the Pell equation 

2 2 i 

x — ny =1 

and a parametric solution of degree d associated to it. 

Proof — Setting 

(2 + V^y = a + bV 3 

and n = 3b 2 , the couple (a, 1) is a fundamental solution of the Pell equation 

2 2 1 

x — ny =1. 

By construction, the integer (U : V (a, b )) is a multiple of d and we can apply theorem 5. □ 

We close the paper by showing that when restricting to square-free integers n, the degree 
of a parametric solution is bounded by 6. This will follow from the elementary result below. 

Lemma 7. — If the integer n is square-free then the group U/ V is cyclic of order dividing 
3. 

Proof. — We have to prove that for any x £ U, we have x 3 £ Z[i Jn], Since for n congruent 
to 2 or 3 modulo 4 we have the identity Ok = Z[y/n], we can assume n = 1 (mod 4), so 
that an element x £ Ok can be uniquely written as x = |(a + b^/n), where a,b £ Z have 
the same parity. Such an element belongs to U if and only if a 2 — rib 2 = 4, in which case we 
find the identities 

8x 3 = a(a 2 + 3 nb 2 ) + b(3a 2 + nb 2 )^/n = 4 a(a 2 — 3) + 46(a 2 — 1 )y/n. 

The integers a and b having the same parity, it follows that a(a 2 — 1) and b{a 2 — 3) are even, 
so that x 3 belongs to Z[y/nJ. □ 

Corollary 8. — Let n be a positive, square-free integer. The degree of a parametric solution 

associated to a fundamental solution of the Pell equation 

2 2 1 

x —ny =1 

is bounded by 6. More precisely, for n congruent to 2 or 3 modulo 4 the degree is 1 or 2 while 
for n congruent to 1 modulo 4 it is equal to 1, 2, 3 or 6. 

Proof. — It is an immediate consequence of theorem 5 combined with lemma 7. □ 
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